o: b 






MICROSCOPIC PREDICTIONS FOR 
CLUSTER-DECAYS 

Mihail MIREA% Aureliu SANDULESCU^'^ Doru Sabin DELION'^''^ 

"'Horia Hulubei National Institute for Physics and Nuclear Engineering, P.O. Box MG-6, 

Bucharest, Romania 



Bioterra university, 81 Garlei str., Bucharest, Romania 
'^Institute for Advanced Studies in Physics, Bucharest, Romania 



Abstract 



The decay dynamical path is determined within the macroscopic-microscopic 
model for the emission of ^'^Ne from ^'^^U. The nuclear shape parametrization 
^ ' is characterized five degrees of freedom. The single particle energies and the 

Ch . nucleon wave functions are obtained within the superasymmetric Woods- 

Saxon two center shell model. It turns out that the cluster decay follows 
»^ , a potential magic valley, starting from the ground state of the parent and 

Q\ ' reaching a configuration of two touching nuclei at scission. A small pocket 

f^ ' in the potential barrier is evidenced, as a result of large shell effects in the 

^ . nascent fragments. The half-life is computed by using several approaches 

h->^ ! for the effective mass. It is shown that the inertia within by the gaussian 

O \ overlap approach gives the closest values to the experimental ones. Half- 

lives for different cluster decays are predicted. The theoretical values are 
compared to various phenomenological estimates. 

• i-H , 

KA , Keywords: Woods-Saxon two center shell model. Cluster decay. Cold 

H ' fission. 



1. Introduction 

The cluster decay was predicted in 80's [jj, I2I, y, 14, la, |6| and experimen- 
tally evidenced in 1984 [a, [ol, llO|, lUl ■ Since then, the spontaneous emission of 
heavy fragments was intensively investigated. An unified phenomenological 
approach of the cluster radioactivity, cold fission and alpha decay as well 
as many body theories were used \T\. Later on, a fine structure in cluster 
emission was also predicted |12i] . Energy spectrum measurements concerning 
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^^C emission from ^^^Ra 13|, llJ] revealed a fine structure with an intense 
branch to the excited state of ^°^Pb. The best agreement between experi- 
ment and theory was obtained by considering the cluster decay as a fission 



process [15|, |16|, treated within the macroscopic-microscopic model. Only the 
experimental hindrance factors were reasonable reproduced. The difference 
between the theoretical and experimental absolute half-lives was subject of 
many orders of magnitude. 

In this work, our aim is to develop a macroscopic-microscopic approach 
to treat in a unitary manner the cluster decay and the fission process, in 
order to reproduce theoretically the half-lives. For this purpose, a fission 
like theory will be used to determine the best sequence of nuclear shapes for 
the cluster decay. In this context, the minimal action principle will be used. 
Two ingredients are needed: the deformation energy of the disintegrating 
system and the nuclear inertia. The half-live of the ^"^Ne emission from ^^^U 
is determined within the WKB approximation and compared to the exper- 
imental values. Such an elaborated study of the fission dynamics in a wide 
range of mass asymmetries could help us to better understand the underly- 
ing physics and to provide an unitary treatment of cluster decay and fission. 
In this respect, the calculations evidenced a pocket shape of the potential 
barrier in the path towards scission. This pocket has a different nature with 
respect to the double potential barrier associated to fission. In fission, the 
second well is obtained as an isomeric state of the parent nucleus. From this 
isomeric state, the single particle levels are strongly rearranged to give the 
asymptotic configuration of two separated fragments. This rearrangement is 
evidenced by the strong positive shell effects of the second barrier. In cluster 
decay the evidenced pocket belongs to a magic valley that leads directly to the 
cluster emission. This valley belongs to a mass asymmetry consistent with 
the formation of the '^^^Pb. Predictions for other decay modes from ^^^U are 
also performed. 

2. Model 

2.1. Fission trajectory 

The calculation addresses ^^Ne cluster emission from ^^^U. The microscopic- 
macroscopic model [17'] is exploited dynamically, by determining the least 
action trajectory. The dynamical analysis of a fissioning nucleus requires 
at least the knowledge of the deformation energy and the effective mass. 
For simplicity, in the macroscopic-microscopic model one assumes that these 



quantities depend upon the shape coordinates. Thus, in our analysis, the 
basic ingredient is the nuclear shape parametrization. The nuclear shape 
parametrization used is given by two ellipsoids of different sizes smoothly 
joined by a third surface obtained by rotating a circle around the symmetry 
axis. Five degrees of freedom characterize this nuclear shape parametrization: 
the elongation, given by the inter-nuclear distance R = Z2 — zi between the 
centers of the ellipsoids, the two deformations of the nascent fragments de- 
noted by their eccentricities e^ = [1 — (ftj/aj)^]"'^/^ (2=1,2), the mass asymmetry 
given by the ratio of major semi-axis r] = 02/01 and the necking parameter 
related to the curvature of the intermediate surface C = s/R^. The quan- 
tity C is used for swollen shapes in the median region, while R3 is used for 
necked shapes. The meaning of the geometric symbols can be understood by 
inspecting Fig. [H A single nucleus and two separated fragments are allowed 
configurations. The determination of the fission trajectory can be obtained 
through a minimization of the action integral in our five- dimensional config- 
uration space, starting with the ground state of the system and ending at 
the exit point of the barrier or in the scission configuration. For fission, such 



calculations were already realized for Th, U [18|, ll9[ and Cf [20 



The probability P of a given channel in a fission process is ruled by an 



exponential factor within the WKB approximation [21 



(1) 

The exponent of the above relation gives the classical action integral of a 
fixed energy along a trajectory in our multidimensional configuration space. 
In the present work, this energy is considered as the ground state energy of 
the parent nucleus. In fission, the trajectory connects the ground state con- 
figuration Ri to the exit point of the barrier Rf. In the cluster decay case, 
the minimization has to be performed up the configuration of two touch- 
ing nuclei, denoted Rf. To this purpose, two ingredients are required: the 
deformation energy V and the tensor of the effective mass B. 

2.2. Deformation energy 

The deformation energy V was obtained by summing the liquid drop 
energy El£,m with the shell and the pairing corrections 6E: 

V = Eldm + 5E . (2) 




Figure 1: Nuclear shape parametrization. 



The macroscopic energy Elom is obtained in the framework of the Yukawa 



plus - exponential model [22| extended for binary systems with different 
charge densities [23| as detailed in Ref. [24 



Eldm = En + Ec + AEc + Ey + Ew + Eao , (3) 

where 

E„ = --^ / /■ (^ - 2) =4Z2) '''M'r, . (4) 

87rV^a4 JvJv^ a / ^ 

is the nuclear term, 

Ec = U [ '-^^i^^^i'nSr, . (6) 

is the Coulomb energy for a constant electric density pe, 



If [1 

2 Jv Jv n 

is a diffuseness correction to the Coulomb potential as described in Ref. 



AEc = -^ — exp-^-/-^ (l + ^) d'nd'r^ , (6) 

2 Jv Jv ^12 ^ 2a / 



22[, and i^y is the volume energy. In the finite range droplet model j25 
there are also several terms that are not deformation dependent. The most 
important are the Wigner term Ew and the A^ energy Eaq- Both were 
added in the potential energy. These terms don't have a shape dependence 
and were neglected in previous investigations. In the following, these terms 
vanish in the overlapping region up to an elongation close to the scission 
configuration Rf. They are supposed to reach linearly their final value for a 
distance larger than 2 fm between the surfaces of the separated fragments, 
where the nuclear forces vanish. In the previous definitions pe are charge 
densities and ri2 =\ fi — r2 \- The charge densities are considered constant 
inside the volume of the nucleus in Eq. (|5]) and zero outside. Another 
degree of freedom has to be introduced here, namely the charge asymmetry. 
Usually [2^1] , one considers that pe is charge density of the parent nucleus for 
an elongation smaller than 0.7 Rt, where Rt is the elongation characterizing 
the configuration of two touching fragments, and varies linearly up to the 
final values of the two nascent nuclei at scission. This charge equilibration 



procedure is described in Ref. [26 



2.3. Shell effects 

The shell effects 5E are obtained as a sum between the shell 5U and pair- 
ing 5P corrections. In this context the Strutinsky procedure was used. These 
corrections represent the varying part of the total binding energy caused by 
the intrinsic structure. In calculating the pairing effect, constant values of 
the pairing matrix elements are computed separately for the parent and the 



two fragments. A renormalization procedure [2l[ in the BCS theory that 
depends on the energy level distribution and a smoothed gap distribution is 
used to obtain the values of the pairing matrix elements associated to each 
fragment issued in the the reaction. To have a smooth transition between 
the shell effects that characterize the parent nucleus to those associated to 
the configuration given by two separated nuclei a simple approximation was 
used. The shell and pairing effects are computed by using the whole level 
scheme from the spherical nuclear shape up to the elongations i?j„t=8.5 fm. 
Beyond Rfin=12 fm the two fragments are separated and there are no ambi- 
guities in determining the wave functions located in the potential well of the 
daughter or in that of the emitted nucleus by using the procedure described 



in Ref. [27|. Thus, beyond Rfin the shell and pairing corrections are de- 
termined separately for each nucleus by using the appropriate level schemes 
for the parent and the daughter nuclei. Between these two configurations 
located at Rint and Rfin, we used a gradual interpolation of the type: 

Uk = Uk + {ul'"- - Uk){R - Rint) /{Rfin " Rint) (7) 

Vk = Vk + (t;f " - Vk){R - Rint) /{Rfin - Rint) , (8) 

where Uk and Vk are the vacancy and occupation amplitudes, respectively, 
calculated for the whole level scheme, while m{™ and f^*"" are their values 
computed separately for each nascent fragment at the elongation Rfin- 

2.4- Effective mass 

The deformation energy is a function of the collective parameters, giving 
the generalized forces acting on the nuclear shape. For a complete description 
of the fission process, it is also necessary to know how the nucleus reacts to 
these generalized forces. This information is contained in the effective mass 



of the system [2l|. The most used approach to calculate the inertia is the 



cranking model. Recently, the cranking model was generalized by taking into 



account the intrinsic excitation produced during the fission process itself |28 . 



In our investigation, three different approximations to describe the inertia are 



used. First of all, for the determination of the fission trajectory, the mass 
parameters are evaluated microscopically within the cranking model. In order 
to evaluate the half-life, two other approximations are used: the Gaussian 
Overlap Approximation (GOA) and the diabatic Cranking formula by using 
Time Dependent Pairing Equations (CTDPE). Relations concerning these 
models are given in the Appendix. 

2.5. Microscopic potential 

To calculate the inertia and the shell effects, we need a microscopic po- 
tential. The microscopic potential has to be consistent with our nuclear 
shape parametrization. The simplest way to define the mean field it to use 
a semi-phenomenological Woods-Saxon potential. In order to take into ac- 
count nuclear deformations going over to separate shapes, a two-center shell 



model with a Woods-Saxon potential was recently developed [16|. The mean 
field potential is defined by the relation: 

Vo{p,z) = \r—T, (9) 



1 + exp 



A(p,^) 



where A{p,z) is the distance between a point {p,z) and the nuclear surface. 
This distance is measured only along the normal direction on the surface 
and it is negative if the point is located in the interior of the nucleus. Vc 
is the depth of the potential, while a is the diffuseness parameter. In our 
work, the depth is K = Voc[l ± k(A^o — ^o)/(A^o + ^o)] with plus sign for 
protons and minus sign for neutrons, Vqc = 51 MeV, a=0.67 fm, k=0.67. 
Here, Aq, Nq and Z^ represent the mass number, the neutron number and 
the charge number of the parent, respectively. This parametrization, usually 
referred as the Blomqvist-Walhlborn one, is adopted because it provides the 
same constant radius tq for the mean field and the pairing field. It ensures 
a consistency of the shapes of the two fields at hyperdeformations, i.e., two 
tangent ellipsoids. The Hamiltonian is obtained by adding the spin-orbit 
and the Coulomb terms to the Woods-Saxon potential. The eigenvalues are 
obtained by diagonalizing the Hamiltonian in the semi-symmetric harmonic 



two center basis ^, |30|. In this work, the maximal principal quantum 
number is Nmax = 14. The two center Woods-Saxon model will be used 
to compute shell and pairing corrections together with inertia. The two 
center shell model represents a valuable instrument to investigate the role of 
individual orbitals for the treatment of a wide variety of superasymmetric 



disintegration processes, pertaining to cluster- and alpha-decays [l5|, l3l|, l32 



or superheavy elements [33|, |34 
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Figure 2: Family of shapes along the minimal action trajectory for cluster decay. The 
elongation R in fni is marked for each shape. 



3. Results and discussion 

3.1. Action integral minimization 

It is not possible to directly minimize the functional (1), due to the large 
computing time determining the values of the potential energies and of the 
effective masses. Thus, in the relevant configuration space, a small number of 
potential energies and inertia are computed and their interpolated values are 
used in the minimization procedure. First of all, a grid of 691 200 deformation 
values was fixed in our five-dimensional configuration space: 20 values of R 
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Figure 3: Upper part: potential energy surface V as function of the elongation R and the 
mass asymmetry rj — axjai- Lower part, contour plot of the potential energy surface. 
The step between two equipotential curves is 2 MeV. The variations of the coordinates 
£1, 62, and C follow the least action path as function of R. The least action trajectory is 
plotted within a thick curve. 



between fm and and the scission point, 8 values of eccentricities e^ between 
and 0.75, 15 values of the ratio ri = ai/a2 in the interval 1 and 3 and 32 values 
for C between 0.11 and 0.20 fm^^, (5 values for R^ between and 4 fm are 
also added). The deformation energy and the elements of the inertia tensor 
were computed in these selected points. In this way, the pertinent region in 
the configuration space including the possible fission trajectories between the 
ground state and the scission configuration Rt was spanned. The quantities of 
interest in this selected region were obtained by interpolating the calculated 
masses and the energies. The trajectory emerges by minimizing numerically 
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Figure 4: The diagonal components of the main effective masses divided by the reduced 
mass ^ of the ^"^Ne emission. Upper paneh Bn/ jjl (dimensionless with step between two 
curves of 0.2). Middle panel: log(i?^/^) (the dimension of i?,,//^ is fm^ and the step 
between two curves is 0.2). Lower panel: \og{Bc I ^j) (the dimension of Be / iJi is fm'' and 
the step between two curves is 0.5). The least action trajectory is plotted by a thick curve. 
Some values of the mass parameter are marked on the plot. 



the action functional, according to Ref. |36| used to describe the fission 
process. 

The family of nuclear shapes for the cluster decay along the least action 
trajectory are plotted in Fig. O The superasymmetric fission trajectory, as 
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Figure 5: Same as Fig. [3] for a representation in the C and R generalized coordinates. 
For C greater that 0.2 fm~^ one uses the rigth scale for _R3 = 1/C. The step between two 
equipotential lines is 2 MeV. 



well as the landscapes of the deformation energies and the main components 
of the tensor of effective masses, are displayed in Figs. El lU E] and El respec- 
tively. In Fig. [31 the potential energy surface is represented as function of 
the mass asymmetry parameter rj = ai/a2 and the elongation R. The depen- 
dencies of all generalized coordinates as function of R follow the variation 
obtained along the minimal action trajectory. The ground state of the parent 
nucleus is located at an elongation R=4.Q fm with a mass asymmetry rj ^ 1. 
The ratio ai/a2 abruptly changes when the nucleus stars to deform, i.e., when 
the elongation increases. A compatible with the final mass asymmetry value 
is very soon obtained. It is clearly evidenced that the nuclear system follows 
a well behaved valley in the potential energy surface up to the scission config- 
uration located at Rt ~ 10 fm. The scission configuration is approximatively 
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Figure 6: Same as Fig. |4] for a representation in the C (or R^) and R generalized 
coordinates. 



described by two tangent spherical fragments. It is interesting to stress that 
in the case of the fission phenomena, the situation is very different. The 
behavior for fission is displayed for ai/a2 ~1.-1.2 fm, where two fragments of 
comparable sized are formed. It can be seen in Fig. 2 that a double barrier 
occurs. The nucleus, initially in the ground state, is fissioning by penetrating 
a first barrier located at -R=6.5 fm and reaches a second well at R=7 fm. 
The situation is completely different for the cluster decay 35|: as mentioned 
the system follows an energy valley in the deformation energy. When the 
elongation is larger than i?^, the valley in the potential energy surface is ex- 
tended in the external region for the system ^°^Pb+^^Ne. Keeping in mind 
that the trajectory is determined by the interplay between the deformation 
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Figure 7: Neutron level diagram for ^^Ne cluster decay from ^32^ •^^^[f]^ respect to the 
elongation R. 

energy and the effective mass, we plotted in Fig. H] the diagonal components 
Brr, Br,r, and Bcc of the inertia tensor as function of 01/02 and C. The 
trajectory has a smooth variations following the contour lines trying to avoid 
large values of the effective masses. A sudden variation of the trajectory is 
not allowed because any rapid variation is translated into an increase of the 
mass. This effect can be understand by appealing to the formula flA.2p . The 
inertia along a path has a strong dependence on the derivatives of the gener- 
alized coordinates. This behavior confirms the aspects discussed in Ref. I 
in the case of the Werner- Wheeler approach. For each mass asymmetry at 
scission, the effective mass Bjm reaches a value close to the reduced mass. 

In Fig. El the deformation energy is displayed in the plane {C,R). The 
ground state s located at i?=4.6 fm and has a necking parameter C=0.075 
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Figure 8: Proton level diagram for ^^Ne cluster decay from ^^"^U with respect to the 
elongation R. 



fin^ . So, in the ground state the shapes are swollen. These swollen shapes 
are preserved up to R=Q fm. From this value, the necking parameter starts 
to vary abruptly and the shapes become very necked producing the rupture 
at Rt ~ 10 fm. At scission, the configuration of nearly two touching nuclei 
is obtained. The variations of the generalized parameters corresponds to an 
increase of the total inertia. The variation of the effective mass is directly 
connected to this evolution of the nuclear system as it can be seen from Fig. 
[6l where the main components of the effective mass are plotted. The values 
of Bcc are large for small or negatives values of C . Thus, in order to have a 
minimal value of the action integral the nucleus has to rapidly escape from 
this region. Notice that the values of the effective mass in the scission region 
are smaller for C — ?■ oo (or R^ =0 fm). Let us mention that Br,T^ decreases 
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Figure 9: Shell effects SU as function of the elongation R: the thick full curve denotes the 
total shell effect 6E, the thin dashed line gives the proton shell effect, the thin dot dashed 
line represents the proton pairing contribution, the thin dotted line is the neutron shell 
effect, while the full thin curve denotes the neutron pairing contribution. 

with mass asymmetry and Bcc vanishes after scission. 

3.2. Level diagrams 

The neutron and proton single-particle diagrams are calculated along the 
minimal action trajectory, from the ground state of the parent nucleus and 
beyond the formation of two separated fragments. These level schemes are 
plotted in Figs. [7] and [8] for neutrons and protons, respectively. In Fig. [71 at 
i? ~ 0, the parent nucleus is considered spherical. For small deformations, the 
system evolves in a way similar to a Nilsson diagram for prolate deformations. 
In the left side of Fig. [71 the orbitals of the parent nucleus (considered 
spherical) are labeled by their spectroscopic notations. The levels of the 
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Figure 10: (a) The maeroscopic-niieroscopic potential barrier V for the ^"^Ne emission as 
a function of the elongation R is plotted by a thick curve. The thin curve represent the 
Coulomb energy ZxZ^e^ jR — Q relative to the experimental Q-value of the process. The 
dashed line represents the barrier given by the liquid drop model, without shell effects, (b) 
The effective mass along the minimal action trajectory calculated within the scmi-adiabatic 
model (full curve), the GOA (dashed curve) and the cranking model (dot-dashed curve). 



emitted fragment can be identified and we plotted them by thick hnes. Both 
fragments are spherical after scission and their levels are bunched in shells. 
The levels of the light fragment are labeled by their spectroscopic notations. 
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Figure 11: (a) Theoretical potential barriers V for the ^^Ne emission (thick curve), for the 
^*Mg emission (dot-dashed curve) and ^^Ne emission (dashed curve) as function of the 
elongation R. (b) The GOA effective mass along the minimal action trajectory for the 
three processes. The same line types are used as in panel (a). 



In the proton diagram of Fig. [HI a decrease of the single particle energies 
after the scission can be observed due to the Coulomb mutual polarization. 
The energy slope for the light fragment is very large. 
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3.3. Half-life of'^'^Ne emission from '^^'^ U 

The shells corrections are displayed in Fig. |9l Notice that beyond R ^ 6 
fm, the total shell effect monotonically decreases up to the asymptotic final 
value. Thus, the nucleus feels the existence of the ^^^Pb magic numbers in 
the overlapping region, well before scission. The macroscopic-microscopic 
potential is plotted in Fig. |3t^a) by a full line. The ground state is located 
at i? ^ 4.6 fm. Due to shells effects, a small pocket occurs in the barrier at 
R p^ 6.5-7 fm. According to Fig. [21 at this elongation the emitted nucleus 
starts its preformation, that is the shape becomes necked in the median 
region of the nuclear shape. Therefore, this pocked occurs as a result of 
the shell effects of nascent fragments. Beyond this point, the increasing 
influence of the liquid drop barrier attenuates the microscopic effects. The 
liquid drop energy is plotted by a dashed line. The Coulomb interaction 
corrected to the experimental Q-value, i.e. the phenomenological quantity 
Vcou = Z1Z2/R — Q is also plotted in order to test the validity of the 
model. The theoretical and the phenomenological values agree well enough 
in the external part of the barrier, the difference being less than 2 MeV. This 
value agree with the r.m.s. deviations of the macroscopic-microscopic model 



25[. In Fig. MJo) the effective mass along the minimal action trajectory 
is plotted for the three mentioned models: cranking, GOA and CTDPE. 
The three inertia exhibit a similar shell structure. The cranking model gives 
the larger values, the GOA the smaller ones, while the CTDPE has always 
intermediate values. The variations of the inertia are mainly generated by 
the slope variation of the generalized coordinates. 

The half-life of ^^Ne emission was estimated by using a semi-empirical 
formula T1/2 = 0.72 x lO^^^P^^ [s]. We used the three models for the in- 
ertia. The values of log(Ti/2[s]) are 35.66, 22.41 and 29.48. They were are 
obtained for cranking model, the GOA and the CTDPE, respectively. The 
experimental value is 21.06. Thus, a reasonable agreement, less than two 
orders of magnitude, is obtained when the inertia is calculated within the 
GOA. 

3.4. Predictions 

We used this model to predict the best candidates for the cluster emission 
process from ^^^U. In order to follow the potential valley, the same trajectory 
was used up to an elongation i? ^ 8 f m and from this point the mass asym- 
metry was gradually changed to reach different final configurations. The 
barriers and the effective masses for the selected reactions are displayed in 



Table 1: Logarithm of the cluster-decay half-lives from ^"^^U. The first column indicates 
the emitted fragment, the second column gives the Q-value, the third column gives the 
experimental values, the fourth column gives the calculated half-lives within GOA effective 
masses, while the fifth column gives the phenomenological values. 
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Fig. [TT] The barriers for ^^Mg and ^^Ne emissions are larger than those 
obtained for ^^Ne. Thus, the magnitude of the shell effects is larger when 
the daughter is the ^osp^ double magic nucleus. The inertia along the least 
action trajectory in plotted in panel (b). The effective masses show similar 
shell structures, asymptotically reaching their appropriate reduced masses. 
Predictions concerning cluster decays from ^^^U are given in table [T] The 
microscopic results are compared to the experimental data |38| and with the 
those given by the analytical superasymmetric fission model |39|, where the 
penetrabilities are obtained by using a phenomenological correction of the 
liquid drop external barrier taking into account the experimental Q-value of 
each process. Both theoretical estimations agree within one or two orders of 
magnitude. 

4. Conclusions 

In this paper, the cluster decay was considered as a superasymmetric fis- 
sion process and it was treated treated within the macroscopic-microscopic 
approximation. The fragmentation potential in the overlap region was ob- 
tained in conjunction with the minimal action principle. It was shown that 
the cluster decay follows a well behaved valley in the potential energy land- 
scape, connected to the formation of the ^°^Pb. A pocket in the potential, 
due to large shell effects of the nascent ^^^Pb, was evidenced. This minimum 
is canceled by a rapid increase of the macroscopic barrier. This behavior is 
very different with respect to the well known double humped fission barrier, 
where the second minimum occurs as a result of the strong shell effects in a 
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isomeric state of the parent nucleus. These results give a better understand- 
ing of the cluster decay phenomenon and represent a contribution towards 
an unitary treatment of fission and heavy ion emission. Usually, the cluster 
decay is treated by calculating a preformation probability, as in the alpha 
decay |7|, |40|. By using GOA inertia, the calculated half-hfe for the ^^Ne 
emission from 232-[j gj^ows a very good agreement with respect to the experi- 
mental data. Microscopic predictions concerning different decay modes were 
also given. 
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Appendix A. Inertia 

Formulas for the elements of the tensor of inertia are given in this Ap- 
pendix. In the cranking approximation, the inertia associated to two gener- 
alized coordinates g^ and qj is 

2 (z/| If |/i)(/i| f^ \iy) 
Mijiqi, g2, -qn) = ^2^ '^(E +E )^''' ^'^"'"'^ ^ ^^^"^'^ ^ ^"^ ^^'^^ 

where v and /x denote the single particle wave functions, E^, u^ and v^, 
are the quasiparticle energy, the vacancy and occupation amplitudes of the 
state z/, respectively, in the BGS approximation, and Pij is a correction that 
depends on the variation of the pairing gap A and of the Fermi energy A as 
function of the deformation coordinates g^. This correction amount up to 10 
% of the total value of the inertia. The inertia B along a trajectory in the 
configuration space spanned by the generalized coordinates g^ (i = l,n) can 
be obtained within the formula 

n n o r^ 

1=1 j=i 

Two other different approximation are tested, namely, the Gaussian Over- 
lap Approximation (GOA) and the diabatic Cranking formula obtained from 
the Time Dependent Pairing Equations (GTDPE). 
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In the generator coordinate method 4l|, |42|, |43|, |4J] the GOA inertia must 
be calculated separately for proton and neutron working spaces: 

1 2 



M„(p) = 2h'' 



T.U,,, Pyi^Pi^viu^v^ + u^v^) 



EujE, + E^)P'^f^P,^AuuV^ + 



(A.3) 



The quantities P^fj. are given by the next formula that depends on a specific 
trajectory in the collective configuration space 



dqi 



P =Vp (a) — 



dR 



(A.4) 



where 42 



dH u^v^ + u^v^ A fdX e^-XdA\ 



dqi ' "■" ' E^, + E^ 



''^2El V% ' A dqj 

(A.5) 

Here, e^ is the single particle energy of the state v. The inertia along a 
trajectory is given by the relation 

(7„ + 7p)'M,Mp 



B 



^iMr, + ^IM, 



(A.6) 



p^"p 



where the index n stands for the neutrons while the index p is for protons. 
Moreover, 7 is computed separately for protons and neutrons 



7 



y^ PufiPf,u{UuVf, + Uf,l 



(A.7) 



!/,/i 



In the case of the CTDPE [28| semi-adiabatic cranking model, the ele- 
ments of the tensor of inertia are: 



M 



^J 






{E,. - Eo) f '7^>' - "''fg-tJ ) < /i 



^^ I z/ >< z/ I 1^ I u > 



dqi 



{Ey^ - Z]^^,/,,. ^7 - -^0 + Yl-y ^7 



(A.8) 

where ti^ = u^Vi, is the paring moment component, p^, = v^ is the occupation 
probability of the level z/ in the seniority zero state. The energies are of the 
seniority two states are: 



E, 



E pT^^ 



u^ — / ^ H'y '-7 
77^ '^.M 



A |2 
G 






v^\2 



+ eu + t. 



(A.9) 
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where the levels e^ and e^ are blocked, the values of p^^ address occupation 
probabilities for the seniority two states and 

T7 = 2p7e, - 2G(p7)^ + "^ ^^' +J"^ ^ ""'' ( i^ - 1 ) (A.IO) 

The energy Eq is the seniority zero state. All quantities without both in- 
dexes v[i address the seniority zero state. The inertia along the trajectory 
is obtained within Rel.( ]A.6|) . The final value of -B is a sum of the quantities 
obtained for protons and neutrons as in the classical cranking model. 

A comparison between the behavior given by these three approximations 
concerning the fission can be found in Ref. j45l] . 
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